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We compute by numerical integration of the Dirac equation the number of quark- 
antiquark pairs produced in the classical color fields of colliding ultrarelativistic nuclei. 
The backreaction of the created pairs on the color fields is not taken into account. While 
the number of qq pairs is parametrically suppressed in the coupling constant, we find 
that in this classical field model it could even be compatible with the thermal ratio to 
the number of gluons. After isotropization one could thus have quark-gluon plasma in 
. chemical equilibrium. 

a: 

' 1. Introduction 

The initial stages of an ultrarelativistic heavy ion collision are believed to be dominated 
by strong classical color fields. There is a twofold interest in calculating the production 
of quark-qntiquark pairs from these classical fields. Firstly, although heavy quark pro- 
duction is in principle calculable perturbatively, it would be interesting to understand 
whether these strong color fields influence the result. Secondly, being able to compute 
both gluon and quark production in the same framework would give insight into the 
chemical equilibration of the system and test the consistency of the assumption of gluon 
dominance. The number of quark pairs present in the early stages of the system has 
observable consequences in the thermal photon and dilepton spectrum. 

In this talk we shall present first results [ Q] of a numerical computation of quark 
antiquark pair production from the classical fields of the McLerran-Venugopalan (MV) 
model. The equivalent calculation, although in another gauge, has been carried out 
analytically to lowest order in the densities of both color sources ("pp"-case) in Ref. [ 
12] and to lowest order in one of the sources ("pA"-case) in Ref. The corresponding 
calculation in the Abelian theory [ IH E] , of interest for the physics of ultraperipheral 
collisions, can be done analytically to all orders in the electrical charge of the nuclei. 
Quark pair production has also been studied in a related "CGC"-approach in Refs. [El EI 
and in a more general background field in Ref. [|H]. 
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Figure 1. Domains of different time de- Figure 2. Dependence on proper time r of the 

pendences. The fermion amplitude is a number of pairs of one flavor per unit rapidity 

sum of two terms: one with interaction dN/ dy for g 2 fi = 2 GeV and for values of quark 

first with the left moving nucleus, then mass marked on the figure. The lowest curve 

the right moving one, and vice versa, corresponds to g 2 fi = 1 GeV. 
A\ ' are pure gauges and Ai,A v is a 
numerically computed color field. 



2. The numerical calculation 

Our calculation of pair production relies on the numerical calculation of the classical 
background color field in which we solve the Dirac equation. 

In the classical field model the background gluon field is obtained from solving the 
Yang- Mills equation of motion with the classical color source J v given by transverse color 
charge distributions of the two nuclei boosted to infinite energy: 

[D^F^] = ,F = 5» + p {1) ( XT )5(x-) + <T-p (2) (x r )<5(x + ). (1) 

In the MV model [ E] the color charges are taken as random variables with a Gaussian 
distribution 

<p a (x T )/(y T )> = <rVW(x T - y T ), (2) 

depending on the coupling g and a phenomenological parameter /x. The combination 
g 2 H is closely related to the saturation scale Q s . Collisions of two ions were first studied 
analytically using this model in Ref. [EI] and the way of numerically solving the equations 
of motion was formulated in Ref. [Illj. 

Our method of solving the Dirac equation is explained in more detail and the numerics 
tested in a 1+1- dimensional toy model in Ref. [ I12j . The domains of spacetime involved in 
the calculation are illustrated in Fig. E One starts in the the infinite past t — > — oo with 
a negative energy plane wave solution ip{x) = e iq x v(q). The Dirac equation can then be 
integrated forward in time analytically to the future light cone (r 2 = 2x + x~ = 0, > 0) 
because the background field in the intermediate region is a pure gauge. This gives an 
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Figure 3. Dependence of the number of 
quark pairs on quark mass at a fixed proper 
time, r = 0.25 fm, and for two values of 
g 2 fi. 



Figure 4. Dependence of the number 
of quark pairs on g 2 \i at a fixed proper 
time, r = 0.25 fm, and for quark mass 
m = 0.3 GeV. 




Figure 5. Transverse momentum spectrum 
of (anti) quarks for g 2 /i = 2 GeV at a fixed 
proper time, r = 0.25 fm, and for different 
quark masses. 



Figure 6. Transverse momentum spectrum 
of (anti) quarks for quark mass m = 0.3 
GeV and for different g 2 \x at a fixed proper 
time, t = 0.25 fm. 



initial condition for numerically solving the Dirac equation for r > using the coordinate 
system r, z, x^. The resulting spinor wavefunction ^(t,z,x.t) is then projected onto 
positive energy states e~ tp ' x u(p) at time r to obtain the amplitude M T . For times larger 
than the formation time of the quark pair r > 1/ ^m 2 + q^ 2 this amplitude can be 
interpreted as the amplitude for producing quark antiquark pairs. The resulting number 
of quark pairs is shown in Fig. El 

The physical parameters of the calculation are g 2 fi characterising the strength of the 
background field, the nuclear radius Ra and the quark mass m. The dependence on g 2 fi 
and m of the number of pairs at r = 0.25 fm is shown in Figs. El and El The transverse 
momentum spectra of the (anti) quarks as a function of is shown for different quark 
masses and saturation scales in Figs. El and [13 
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3. Discussion 

According to conventional wisdom the initial state of a heavy ion collision is dominated 
by gluons. Assuming that the subsequent evolution of the system conserves entropy 
this would mean ~ 1000 gluons in a unit of rapidity. In the classical field model this 
corresponds to g 2 [i ~ 2 GeV. Our results seems to point to a rather large number 
of quark pairs present already in the initial state. One could envisage a scenario where 
for g 2 ji « 1.3 GeV these 1000 particles could consist of > 400 gluons, > 300 quarks and 
> 300 antiquarks (take the lowest curve from Fig. |2] and multiply by Nf — 3). This would 
be close to the thermal ratio of N g /N q = 64/(21iV f ). 

4. Conclusions 

We have calculated quark pair production from classical background field of McLerran- 
Venugopalan model by solving the 3+1-dimensional Dirac equation numerically in this 
classical background field. We find that number of quarks produced is large, pointing to 
a possible fast chemical equilibration of the system. The mass dependence of our result is 
surprisingly weak and we are not yet able to make any conclusions on heavy quarks until 
studying the numerical issues involved. 
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